The Weibull - log Weibull distribution for interoccurrence times of earthquakes 
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By analyzing the Japan Meteorological Agency (JMA) seismic catalog for different tectonic set- 
. . . tings, we have found that the probability distributions of time intervals between successive earth- 

00 ' 

, quakes — interoccurrence times — can be described by the superposition of the Weibull distribution 

, and the log- Weibull distribution. In particular, the distribution of large earthquakes obeys the 

^ , Weibull distribution with the exponent ai < 1, indicating the fact that the sequence of large 

1^ , earthquakes is not a Poisson process. It is found that the ratio of the Weibull distribution to the 

\ probability distribution of the interoccurrence time gradually increases with increase in the thresh- 

(N ■ old of magnitude. Our results infer that Weibull statistics and log- Weibull statistics coexist in the 

interoccurrence time statistics, and that the change of the distribution is considered as the change 
' fH ' I of the dominant distribution. In this case, the dominant distribution changes from the log- Weibull 

distribution to the Weibull distribution, allowing us to reinforce the view that the interoccurrence 
O ' time exhibits the transition from the Weibull regime to the log- Weibull regime. 



PACS numbers: 91.30.Dk, 91.30.Px, 05.65.-|-b, 05.45.Tp 



I. INTRODUCTION 



' Earthquakes arc phenomena exhibiting great complexity eharaeterized by many empirieal statistieal laws [l|. The 

lO time intervals between successive earthquakes can be classified into two types: interoccurrence times and recurrence 

oo : . n 

, times [2[ . Interoccurrence times are the interval times between earthquakes on all faults in a region, whereas recurrence 

, times are the time intervals between earthquakes in a single fault or fault segment. For a seismologist, recurrence 

t — • . 

■ times denote the interval times of characteristic earthquakes that occur quasi-pcriodically in a single fault. 

' Recently, a unified scaling law of interoccurrence times was reported using the Southern California Q| and worldwide 

^ ' earthquake catalogs where the interoccurrence times were analyzed for the events with the magnitude m above a 

I certain threshold rUc under the following two conditions: (a) earthquakes can be considered as a point process in space 

^ ] and time; (b) there is no distinction between foreshocks, mainshocks, and aftershocks. It has been demonstrated that 
the probability distribution of the interoccurrence time is well-fitted by the generalized gamma distribution. This 
scaling law is obtained by analyzingthe aftershock data Q| and is derived approximately from a theoretical framework 



proposed by Saichev and Sornctte [6[ . Abe and Suzuki showed that the survivor function of the interoccurrence time 
can be described by a power law \A. It has been reported that the sequence of aftershocks and successive independent 
earthquakes is a Poisson process [3, S] ■ Recent works on interoccurrence time statistics are focused on the effect of 
"long-term memory" 
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121 as well as on the determination of the distribution function. However, the effect of 
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changing the threshold of magnitude on the interoccurrence time statistics has not been discussed fully in the existing 
literature. 

In this work, we trace the change in interoccurrence time statistics produced by varying of the cutoff magnitude 
mc, which has not been studied previously. This study aims to infer the interoccurrence time statistics for middle or 
big mainshocks. It demonstrates that the interoccurrence time distribution is described by the superposition of the 
Weibull distribution and the log-WeibuU distribution. In particular, the distribution of large earthquakes follows the 
Weibull distribution with exponent ai < 1, indicating that the sequence of large earthquakes is not a Poisson process. 
In addition, the ratio of the Weibull distribution to the interoccurrence time distribution gradually increases as riic 
is increased. Our results lead to the conclusion that the interoccurrence time statistics follow Weibull statistics and 
log- Weibull statistics closely, and the change of the distribution can be interpreted as the change of a predominant 
distribution, i.e., the predominant distribution changes from the log- Weibull distribution to the Weibull distribution 
when rric is increased, so that the interoccurrence time statistics exhibit the Weibull - log Weibull transition. 



II. DATA ANALYSIS AND APPLICABLE DISTRIBUTIONS 



In this study, we investigate interoccurrence time statistics using the earthquake catalog made available by the Japan 
Meteorological Agency (JMA) [l3|. This catalog contains occurrence times and hypocenter locations of earthquakes 
having a magnitude m greater than 2.0, and covers the region spanning 25°-50° N latitude and 125°-150° E longitude 
(see Fig. [T](a)). In this study, we use the data from January 1, 2001 to October 31, 2007. 



Quary i, 
lB[l2 



Our method is similar to the that of previous works 

UMM (see Fig. [ID; 

1. We divided the spatial areas into a window of L degrees in longitude and L degrees in latitude. 

2. For each bin, earthquakes with magnitude m above a certain cutoff magnitude uic were considered. 

3. We analyzed the interoccurrence times and then performed the data fitting in the time domain, t > h (day). 

For each bin. we analyzed interoccurrence times using at least 100 events to avoid statistical errors, h and L are taken 
to be 0.5 and 5, respectively. As shown in Fig. [T] (a), we investigated the interoccurrence time statistics for different 

nn 

14 regions. Aftershocks might be excluded from the study based on the information from previous studies [J, |9|. 

One of our main goals in this study is to determine the distribution hmction of the interoccurrence time. Here, 
we will focus our attention on the applicability of the Weibull distribution P^, the log- Weibull distribution Pi^ [2ll |. 
the power law P^ow Qi the gamma distribution Pgam (in the case of (5 = 1 in the paper [J]), and the log normal 
distribution Pin [22], which are defined as 

ai — 1 
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where a^, (3i, and h are constants and characterize the distribution. T{x) is the gamma function, i stands for an index 
number; i = 1,2,3,4, and 5 correspond to the Weibull distribution, the log- Weibull distribution, the power law, the 
gamma distribution, and the log normal distribution, respectively. 
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FIG. 1: (a): A map around Japan where we have carried out analysis the JMA catalog. The number of each bin stands for the 
number of earthquakes from 01/01/2001 to 10/31/2007. (b) Illustration of interoccurrence times for three different threshold 
of magnitude, mci,mc2, and mcs. In our work, we focus on interoccurrence times in the time domain, t > h, corresponding to 
o in the figure. 



The WcibuU distribution is well known as a description of the probability distribution of failure-occurrence times. 



mechanics of rocks 
the ID Q and 2D 
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In seismology, the distribution of ultimate strain l^l, the recurrence time distribution [isl . [l^, and the damage 



show the WeibuU distribution. In numerical studies, the recurrence time distribution in 



18j spring-block model, and in the "Virtual California model" [19[ also exhibit the WeibuU 
distribution. It is known that for ai = I and ai < 1, the tail of the WeibuU distribution is equivalent to the 
exponential distribution and the stretched exponential distribution, respectively. The log- WeibuU distribution is 
constructed by a logarithmic modification of the cumulative distribution of the WeibuU distribution. In general, the 
tail of the log- WeibuU distribution is much longer than that of the WeibuU distribution. As for Q!2 = 1 the log- WeibuU 
distribution is equal to the power law. It has been reported that the log- WeibuU distribution can be derived from the 
chain-reaction model proposed by Huillet and Raynaud [21 1 . 

In order to detect whether a specific distribution is preferred, we used three different goodness-of-fit tests: the root 
mean square (rms) test, the Kolmogorov-Smirnov (KS) test, and the Anderson-Darling (AD) test. Firstly, we use the 
rms test. The rms value is defined as 



where Xi are actual data of the distribution and x'^ are predicted data obtained from the ideal curve, n and k stand 
for the number of data point and fitting parameters, respectively. In this study, the rms value is calculated using 
the cumulative distribution function (cdf) and the WeibuU plot. As readily known, the most appropriate distribution 
shows the smallest rms value. Secondly, the KS test is performed. The Kolmogorov-Smirnov statistic Dks is defined 
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FIG. 2: Distribution of interoccurrence time at Okinawa region for different rric and distribution functions, (a) rric — 4.5, (b) 
TTic = 3.0, and (c) tUc = 2.0. Inset figure of (a) shows the Weibull plot and (b) and (c) represent the log-log scale of the figures. 

TABLE I: The results of the rms value, Dks, A^, and optimal parameters for different distribution functions in Fig. [2] (a). The 
error bars are 95 per cent confidence limits. 





distribution 






rms (cdf) 


Dks 




rms (WeibuU-plot) 




P» {i = 1) 


0.82± 0.01 


17.5± 0.28 


0.014 


0.039 


0.89 


0.065 


nic = 4.5 


P» (i = 2) 


2.95± 0.12 


32.4 ± 1.08 


0.028 


0.107 


1.90 


0.099 


127 events 


Ppow — 3) 


1.46 ± 0.04 


0.99± 0.24 


0.112 


0.384 


14.8 


0.19 


91 data points 


Pgam (« = 4) 


0.94± 0.01 


17.6 ± 0.44 


0.022 


0.066 


1.51 


0.098 




Pin {i = 5) 


2.37± 0.04 


1.22 ± 0.05 


0.027 


0.088 


1.06 


0.093 



by 

Dks = T^a.x\yi-y'i\, (3) 

z 

where yi and y[ mean the actual data of the cumulative distribution function and the data estimated from the fitted 
distribution, respectively. It is well recognized that the preferred distribution has the smaUest value of Dks- Finally, 
the AD test is used. This test gives more weight to the tails of the distribution than the KS test. The Anderson-Darling 
statistic is defined as 

1 ^ 

A' = -N-- ^(2* - l)pn(F(zO) + ln(l - F(z„_.+i))], (4) 
1=1 

where N is the sample number, F{x) is the cumulative distribution function, and Zi is the ith interoccurrence time. 
Note that z is put in order. The most suitable distribution exhibits the smallest A^. We applied these goodness-of-fit 
tests by use of the five distributions which were mentioned before. 

III. RESULTS 

The cumulative distributions of the interoccurrence time for different nic in Okinawa region (125°-130°E and 
25°-30°N) and in Chuetsu region (135°-140°E and 35°-40°N) are displayed in Figure [2] and El respectively. The total 
number of earthquakes were 12024 in Okinawa and 13678 in Chuetsu. 
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FIG. 3: Cumulative distribution of interoccurrence time for Chuetsu area at different rric and distribution function, (a), (b), 
and (c) represent interoccurrence time when rric — 4.0, rric — 3.0, and rric — 2.0, respectively. Inset figure of (a) shows the 
WeibuU plot, and (b) and (c) represent log-log scale of the figures. 

TABLE II: The results of the rms value, Dks, A^, and optimal parameters for different distribution functions in Fig. [3] (a). 
The error bars are 95 per cent confidence limits. 





distribution 




ft 


rms (cdf) 


Dks 




rms (WeibuU-plot) 




P» (i = 1) 


0.79± 0.02 


9.35 ± 0.28 


0.019 


0.034 


3.17 


0.083 


rric = 4.0 


Ply, {i = 2) 


2.36 ± 0.18 


16.9 ± 1.08 


0.038 


0.106 


3.71 


0.134 


231 events 


Ppow — 3) 


1.52 ± 0.06 


0.68 ± 0.14 


0.100 


0.251 


14.1 


0.227 


57 data points 


Pgam (i = 4) 


0.95 ± 0.02 


9.70 ± 0.48 


0.034 


0.072 


3.41 


0.077 




Pin {i = 5) 


1.71 ± 0.06 


1.27 ± 0.07 


0.033 


0.056 


1.33 


0.124 



At first, we focus on the interoccurrence time distribution for large TOc- We have tried four statistical tests, the 
rms (cdf), the KS test, the AD test, and the rms (WeibuU-plot) test, whose results are shown in TableUfor Okinawa 
and Table [Til for Chuetsu. For Okinawa, we can certify that the most suitable distribution is the Weibull distribution 
in all tests. In general, there is a possibility that the preferred distribution is not unique but depends on the test we 
use. However, the results obtained in Table [T] seems to support that the Weibull distribution is the most appropriate 
distribution in this case. As for Chuetsu, by four tests, the preferred distribution is suited to be the Weibull distribution 
as shown in Table IHl where the Weibull distribution is the most prominent distribution in the two tests (rms (cdf) 
and KS test), although the Weibull distribution is not the most appropriate distribution in the remaining two tests 
(AD test and rms (WeibuU-plot)). Thus, we reinforce the view that the Weibull distribution is preferred. Hereinafter 
the preferred distribution function is evaluated by use of the rms (cdf) test. 

However, the fitting accuracy of the Weibull distribution becomes worse with a gradual decrease in rric- We 
now propose a possible explanation that states that "the interoccurrence time distribution can be described by the 
superposition of the Weibull distribution and another distribution, hereafter referred to as the distribution X, Px", 

P{t)=pxP^ + (1-p)xPx, (5) 

where p is a parameter in the range, < p < 1 and stands for the ratio of divided by P{t). The interoccurrence 
time distribution obeys the Weibull distribution for p = 1, whereas it follows the distribution X for p = 0. In 
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TABLE III: The interoccurrence time statistics of earthquakes in Okinawa region. The error bars are 95 per cent confidence 
limits. 



rric Distribution X Weibull distribution Distribution X Weibull rate rms-value 



Region 


index i 






P 


[xlO ^] in rms 


4.5 

Okinawa 


Pin, {i = 2) 
Ppow {i — 3) 

Pgam {i = 4) 

Pin {i = 5) 


0.82 ± 0.01 17.5 ± 0.28 
0.82 ± 0.01 17.5 ± 0.28 
0.82 ± 0.01 17.5 ± 0.28 
0.82 ± 0.01 17.5 ± 0.28 


- 


1 
1 
1 
1 


14 -4.27 
14 -4.27 
14 -4.27 
14 -4.27 


4.0 

Okinawa 


Pin, {i = 2) 

-Ppow ('^ — 3) 
Pgam (l = 4) 

Pin {i = 5) 


0.91 ± 0.01 8.28 ± 0.12 
0.91 ± 0.01 8.28 ± 0.12 
0.91 ± 0.01 8.28 ± 0.12 
0.91 ± 0.01 8.28 ± 0.12 




1 
1 
1 
1 


11 -4.51 
11 -4.51 
11 -4.51 
11 -4.51 


3.5 

Okinawa 


Ph. = 2) " 

Ppow — 3) 
Pgam (i = 4) 
Pin (i = 5) ' 


1.09 ± 0.03 3.74 ± 0.13 
1.09 ± 0.03 3.74 ± 0.14 
1.09 ± 0.03 3.74 ± 0.14 
1.09 ± 0.03 3.74 ± 0.14 


1.87 ± 0.25 4.73 ± 0.70 
1.81 ± 0.07 0.64 ± 0.07 

0.86 ± 0.01 0.93 ± 0.01 


0.78 ± 0.03 
0.95 ± 0.02 
1 

0.66 ± 0.06 


5.1 -5.28 
7.3 -4.92 
12 -4.42 
5.7 -5.17 


3.0 

Okinawa 


Pi» (j = 2) " 

Ppow — 3) 
Pgam = 4) 
Pin (i = 5) ' 


1.43 ± 0.04 1.81 ± 0.04 
1.38 ± 0.04 1.64 ± 0.02 
1.38 ± 0.04 1.64 ± 0.02 
1.38 ± 0.04 1.64 ± 0.02 


1.34 ± 0.06 2.46 ± 0.08 
2.20 ± 0.08 0.55 ± 0.02 

0.19 ± 0.006 1.46 ± 0.01 


0.60 ± 0.02 
0.79 ± 0.02 
1 

0.08 ± 0.08 


3.8 -5.57 
8.8 -4.73 
17 -4.07 
6.4 -5.05 


2.5 

Okinawa 


P™ (i = 2) " 

Ppow — 3) 
Pgam (l = 4) 

Pin {i = 5) 


1.76 ± 0.04 1.16 ± 0.02 
1.90 ± 0.09 1.02 ± 0.01 
1.90 ± 0.05 1.02 ± 0.008 


1.23 ± 0.03 1.67 ± 0.02 
2.83 ± 0.08 0.51 ± 0.01 
1.09 ± 0.03 0.94 ± 0.04 
-0.19 ± 0.008 0.53 ± 0.01 


0.47 ± 0.01 
0.59 ± 0.04 
0.99 ± 0.09 



3.0 -5.81 
9.8 -4.63 
23 -3.77 
12 -4.42 


2.0 

Okinawa 


P^. (j = 2) 

Ppow — 3) 
Pgam (« = 4) 
Pn (i = 5) 


1.75 ± 0.06 0.78 ± 0.02 
2.56 ± 0.18 0.77 ± 0.01 
2.56 ± 0.18 0.77 ± 0.01 


1.18 ± 0.03 1.43 ± 0.02 
3.60 ± 0.09 0.48 ± 0.004 
1.09 ± 0.06 0.68 ± 0.05 
-0.41 ± 0.01 0.39 ± 0.03 


0.42 ± 0.02 
0.39 ± 0.04 
0.96 ± 0.08 



2.3 -6.07 
7.3 -4.92 
25 -3.09 
15 -4.20 



"We used the parameter estimation procedure (A). 
''We used the parameter estimation procedure (B). 

this study, the log- Weibull distribution, the power law, the gamma distribution, and the log normal distribution are 
candidates for the distribution X. 

We shall explain the parameter estimation procedure; (A) the optimal parameters are estimated so as to minimize 
the differences between the data and the test function by varying five parameters, ai, a^, and p. If there is 
a parameter, where C^,, the ratio of the standard deviation divided by the mean exceeds 0.1, another estimation 
procedure is performed. (B) the Weibull parameters, ai and /3i, and the distribution X parameters, and /3i, are 
optimized dependently and then p is estimated. 

On the basis of this hypothesis, the fitting results of P{t) are listed in Table Hill for Okinawa region and in Table HVl 
for Chuetsu region. We can assume that the Weibull distribution is the fundamental distribution, because p = 1 
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TABLE IV: The interoccurrence time statistics of earthquakes in Chuetsu area. The error bars are 95 per cent confidence 
hmits. 



rUc Distribution X WeibuU distribution Distribution X Weibull rate rms-value 



Region 


index i 


ai /3i 


ai Pi 


P 


[xlO ^] In rms 


4.0 

Chuetsu 


Pin, {i = 2) 

Ppow — 3) 
Pgam {i = 4) 

Pin {i = 5) 


0.79 ± 0.02 9.35 ± 0.28 
0.79 ± 0.02 9.35 ± 0.28 
0.79 ± 0.02 9.35 ± 0.28 
0.79 ± 0.02 9.35 ± 0.28 




1 
1 
1 
1 


19 -3.96 
19 -3.96 
19 -3.96 
19 -3.96 


3.5 

Chuetsu 


Pin, (i = 2) " 

Ppow — 3) 
Pgam (i = 4) 
Pn = 5) " 


0.85 ± 0.007 4.56 ± 0.03 
0.85 ± 0.007 4.56 ± 0.03 
0.85 ± 0.007 4.56 ± 0.03 
0.85 ± 0.007 4.56 ± 0.03 


1.96 ± 0.03 8.16 ± 0.19 
1.66 ± 0.03 0.60 ± 0.05 
0.96 ± 0.004 4.61 ± 0.05 
1.05 ± 0.02 2.30 ± 0.04 


0.86 ± 0.02 
0.96 ± 0.01 
0.92 ± 0.07 
0.78 ± 0.04 


6.2 -5.08 
7.1 -4.95 
7.6 -4.88 
6.4 -5.05 


3.0 

Chuetsu 


Pn, (i = 2) " 
-Ppow — 3) 
Pgam (i = 4) 

Pn (i = 5) " 


1.08 ± 0.02 2.17 ± 0.02 
1.08 ± 0.02 2.17 ± 0.02 
1.08 ± 0.02 2.17 ± 0.02 
1.08 ± 0.02 2.17 ± 0.02 


1.99 ± 0.16 5.35 ± 0.39 
1.97 ± 0.04 0.53 ± 0.02 

0.40 ± 0.02 0.92 ± 0.04 


0.82 ± 0.04 
0.93 ± 0.009 
1 

0.65 ± 0.06 


3.9 -5.55 
5.0 -5.30 
6.5 -5.04 
3.7 -5.60 


2.5 

Chuetsu 


Pn, (i = 2) " 
Ppow ('^ — 3) 
Pgam (i = 4) 

Pn {i = 5) 


1.47 ± 0.03 1.24 ± 0.02 
1.55 ± 0.06 1.17 ± 0.02 
1.55 ± 0.06 1.17 ± 0.02 


1.20 ± 0.04 1.91 ± 0.04 
2.52 ± 0.08 0.50 ± 0.02 
1.03 ± 0.02 1.10 ± 0.05 
-0.10 ± 0.006 0.64 ± 0.008 


0.59 ± 0.01 
0.69 ± 0.03 
0.99 ± 0.09 



2.4 -6.03 
6.8 -4.99 
16 -4.14 
5.1 -5.28 


2.0 

Chuetsu 


Pn, (i = 2) " 
Ppow — 3) 
Pgam (i = 4) 

Pn {i = 5) 


1.77 ± 0.08 0.78 ± 0.02 
2.43 ± 0.18 0.79 ± 0.01 
2.43 ± 0.18 0.79 ± 0.01 


1.20 ± 0.04 1.47 ± 0.02 
3.47 ± 0.12 0.48 ± 0.008 
1.19 ± 0.11 0.69 ± 0.04 
-0.40 ± 0.005 0.82 ± 0.02 


0.47 ± 0.02 
0.47 ± 0.05 
0.97 ± 0.11 



2.8 -5.88 
8.0 -4.83 
26 -3.65 
14 -4.26 



"We used the parameter estimation procedure (A). 
''We used the parameter estimation procedure (B). 

for large rric, which means that the effect of the distribution X is neghgible. As observed in Table Hill and HVl the 
log- Weibull distribution is the most suitable distribution for the distribution X. Thus, we find that the interoccurrence 
times distribution can be described by the superposition of the Weibull distribution and the log- Weibull distribution, 
namely, 

P(t) =px P^ + (l-p) X (6) 

P(r) is controlled by five parameters, ai, a2, /?2, and p. 

The result of fitting parameters of P{t) for different regions are listed in Table fVl It was found that the interoccur- 
rence time distribution of earthquakes with large iric obeys the Weibull distribution with the exponent ai < 1, which 
has been observed previously. As shown in Tables IIIIl IIVI and |Vl we stress the point that the distribution function 
of the interoccurrence time changes with varying TOc. This indicates that the interoccurrence time statistics basically 
contain both Weibull and log- Weibull statistics, and the change of distribution function can be interpreted as the 
change of a dominant distribution. In this case, the dominant distribution of the interoccurrence time changes from 
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FIG. 4: The crossover magnitude m** map around Japan, m** depends on the region and ranges from 2.4 to 4.3. 



the log-WeibuU distribution to the WeibuU distribution with rric increased. Thus, the point that the interoccurrence 
time statistics exhibit transition from the WeibuU regime to the log- WeibuU regime is reinforced. It is noted that a 
crossover magnitude from the superposition regime to the WeibuU regime, denoted by m**, depends on the spatial 
area on which we have focused in this study. We display m** map around Japan in Fig. |4l m** ranges from 2.4 
(125°E-130°E and 30°N -35°N) to 4.3 (130°E-135°E and 30°N-35°N, and 145°E-150°E and 40°N-45°N). Comparing 
Fig. [T] (a) with Fig. [H we found that a WeibuU - log WeibuU transition occurs in all region where we conducted. 



IV. DISCUSSION 



In our study of the size-dependence of the interoccurrence time statistics, the window size L varied from 3° to 25°. 
We used the data covering the region 140°-143° E and 35°-38° N for L = 3, 140°-145° E and 35°-40° N for L = 5, 
140°-150° E and 35°-45° N for L = 10, and 125°-150° E and 25°-50° N for L = 25. For L ^ 25, the data covers the 
entire region of the JMA catalog. The result of fitting parameters of P{t), the crossover magnitude m**, and the rms 
value are listed in Table IVll We have demonstrated that the WeibuU exponent ai is less than unity and the WeibuU 
- log WeibuU transition occurs in all cases, m** depends on L, namely m** = 3.9 for L = 3, rn** = 4.0 for L = 5, 
TO** — 4.2 for L — 10, and m** ~ 5.0 for L ~ 25. Therefore we can conclude that the interoccurrence time statistics 
presented are valid from L = 3toL = 25. 

Finally, we compared our results with those of the previous studies. The unified scaling law shows a generalized 
gamma distribution which is approximately the gamma distribution, because S in Corral's paper Q is close to unity 
(d ~ 0.98 ± 0.05). For a long time domain, this distribution decays exponentially, supporting the view that an 
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TABLE V: Summary of the interoccurrence time statistics for different regions in Japan. The error bars are 95 per cent 
confidence limits. 

Weibull distribution log-Weibull distribution Weibull rate 



Region 




Ql 


/3i [day] 


a2 


P2 [day] 


P 




2.0 " 


1.26 ± 0.02 


1.59 ± 0.05 


1.43 ± 0.08 


2.07 ± 0.05 


0.59 ± 0.03 


(A) 


2.5 " 


1.31 ± 0.02 


1.90 ± 0.02 


1.45 ± 0.09 


2.09 ± 0.04 


0.67 ± 0.02 


140°E-145°E 


3.0 ' 


1.13 ± 0.02 


1.80 ± 0.02 


1.51 ± 0.16 


3.20 ± 0.10 


0.92 ± 0.04 


25°N-30°N 


3.5 ' 


1.11 ± 0.01 


2.49 ± 0.01 


1.80 ± 0.13 


1.47 ± 0.01 


0.95 ± 0.02 




4.0 


0.97 ± 0.01 


4.45 ± 0.02 


- 


- 


1 




4.5 


0.89 ± 0.02 


11.3 ± 0.19 


- 


- 


1 




2.0 ' 


4.34 ± 0.73 


0.58 ± 0.02 


1.04 ± 0.08 


1.21 ± 0.01 


0.40 ± 0.06 




2.5 " 


2.91 ± 0.28 


0.79 ± 0.02 


1.10 ± 0.05 


1.43 ± 0.02 


0.31 ± 0.04 


(B) 


3.0 " 


1.54 ± 0.04 


1.14 ± 0.04 


1.37 ± 0.06 


1.76 ± 0.06 


0.62 ± 0.02 


140°E-145°E 


3.5 ' 


1.22 ± 0.02 


1.82 ± 0.02 


1.59 ± 0.14 


3.27 ± 0.02 


0.85 ± 0.04 


35°N-40°N 


4.0 


0.94 ± 0.009 


3.45 ± 0.02 


- 


- 


1 




4.5 


0.82 ± 0.02 


7.69 ± 0.14 


- 


- 


1 




5.0 


0.84 ± 0.03 


18.7 ± 0.55 






1 




2.0 " 


2.17 ± 0.04 


0.76 ± 0.006 


1.11 ± 0.02 


1.39 ± 0.02 


0.52 ± 0.02 


(C) 


2.5 " 


1.82 ± 0.05 


1.04 ± 0.02 


1.22 ± 0.04 


1.66 ± 0.02 


0.50 ± 0.02 


135°E-140°E 


3.0 " 


1.37 ± 0.03 


1.71 ± 0.03 


1.33 ± 0.05 


2.24 ± 0.07 


0.60 ± 0.02 


30°N-35°N 


3.5 ' 


1.01 ± 0.01 


3.63 ± 0.06 


1.90 ± 0.18 


4.91 ± 0.55 


0.87 ± 0.02 




4.0 ' 


0.88 ± 0.01 


8.24 ± 0.09 


2.52 ± 0.10 


15.1 ± 0.47 


0.95 ± 0.08 




4.5 


0.93 ± 0.01 


21.6 ± 0.19 






1 



"We used the parameter estimation procedure (A). 
''We used the parameter estimation procedure (B). 

TABLE VI: The interoccurrence time statistics for different system size L. The error bars are 95 per cent confidence limits. 



L 


Region 


m** 


m-c 


ai 


/3i [day] 


rms 


In rms 


L = 3 


140°-143° E and 35°-38° N 


3.9 


4.6 


0.88 ± 0.02 


19.4 ± 0.36 


0.011 


-4.51 


L = 5 


140°-145° E and 35°-40° N 


4.0 


4.7 


0.75 ± 0.03 


10 ± 0.38 


0.014 


-4.27 


L = 10 


140° -150° E and 35°-45° N 


4.2 


4.9 


0.94 ± 0.01 


8.36 ± 0.08 


0.0077 


-4.87 


L = 25 


125°-150° E and 25°-50° N 


5.0 


5.7 


0.93 ± 0.03 


17.8 ± 0.40 


0.021 


-3.86 



earthquake is a Poisson process. However, we have demonstrated that the Weibull distribution is more appropriate 
than the gamma distribution because the rms value obtained from the Weibull distribution is smaller than that from 
the gamma distribution. In addition, for large rric, the probability distribution in a long time domain is similar to 
the stretched exponential distribution because ai is less than unity, suggesting that earthquakes obey the long-tail 
distribution. We stress the point that the probability distribution changes by varying TOc, supporting the view that a 
transition occurs from the Weibull regime to the log-Weibull regime, which has not been reported previously. Recently, 
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2000 4000 6000 8000 500 1000 1500 20 40 60 80 100 



Interoccurrence time Interoccurrence time Interoccurrence time 



FIG. 5; The interoccurrence time distribution for different rric obtained from the spring-block modeL More details are studied in 
Ref. The WeibuU regime, the superposition regime, and the log-WeibuU regime are shown in (a), (b), and (c), respectively. 
Inset figure of (a) and (c) are the WeibuU-plot and log-WeibuU-plot, respectively. 



Abaimov et al. showed that the recurrence time distribution is also well-fitted by the Weibull distribution 31 rather 

. . . n 

than the Brownian passage time (BPT) distribution [22|] and the log normal distribution. Taken together, we infer 
that both the recurrence time statistics and the interoccurrence time statistics show the Weibull distribution. 

In this study, we propose a new insight into the interoccurrence time statistics, stating that the interoccurrence 
statistics exhibit the Weibull - log Weibull transition. This stresses that the distribution function can be described by 
the superposition of the Weibull distribution and the log- Weibull distribution, and that the predominant distribution 
function changes from the log- Weibull distribution to the Weibull distribution as rric is increased. Note that there is 
a possibility that a more suitable distribution might be found instead of the log- Weibull distribution. However, since 
our results are also obtained by analyzing the Southern California and the Taiwan earthquake catalogs as well [2^ 
we believe that the log-WeibuU distribution is the best. Furthermore, the Weibull - log Weibull transition can be 
extracted more clearly by analyzing synthetic catalogs produced by the spring-block model (Fig. [5]) [fsj . 



V. CONCLUSION 



In conclusion, we have proposed a new feature of interoccurrence time statistics by analyzing the JMA earthquake 
catalogs for different tectonic conditions. We found that the probability distribution of the interoccurrence time can 
be described clearly by the superposition of the Weibull distribution and the log- Weibull distribution. Especially for 
large earthquakes, the interoccurrence time distribution obeys the Weibull distribution with the exponent ai < 1, 
indicating that a large earthquake is not a Poisson process but a phenomenon exhibiting a long-tail distribution. As 
the threshold of magnitude rric increases, the ratio of the Weibull distribution to the interoccurrence time distribution 
p gradually increases. Our findings support the view that the Weibull statistics and log- Weibull statistics coexist in 
the interoccurrence time statistics. We interpret the change of distribution function as the change of the predominant 
distribution function; the predominant distribution changes from the log- Weibull distribution and the Weibull distri- 
bution when rric is increased. Therefore, it is concluded that the interoccurrence time statistics exhibit the Weibull - 
log Weibull transition. We believe that this work is a first step toward a theoretical and geophysical understanding 
of this transition. 
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